Effects of anisotropy on thermal entanglement 
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We study the thermal entanglement in the two-qubit anisotropic XXZ model and the Heisen- 
berg model with Dzyaloshinski-Moriya (DM) interactions. The DM interaction is another kind of 
anisotropic antisymmetric exchange interaction. The effects of these two kinds of anisotropics on 
the thermal entanglement are studied in detail for both the antiferromagnetic and ferromagnetic 
cases. 
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Recently the concept of thermal entanglement was in- 
troduced and studied within one-dimensional isotropic 
Heisenberg model [Q. The state of the system de- 
scribed by the Hamiltonian H at thermal equilibrium is 
p{T) = exp{-^) /Z, where Z =Tr[exp(-^)] is the 
partition function and k is the Boltzmann's constant. As 
p(T) represents a thermal state, the entanglement in the 
state is called the thermal entanglement |l| . 

For two-qubit isotropic Heisenberg model there exists 
thermal entanglement for the antiferromagnetic case and 
no thermal entanglement for the ferromagnetic case [Q. 
While for the XY model the thermal entanglement ap- 
pears for both the antiferromagnetic and ferromagnetic 
cases . It is known that the isotropic Heisenberg model 
and the XY model are special cases of the anisotropic 
Heisenberg model (see Eq.(^)). So it is worth to study 
the thermal entanglement in the anisotropic models and 
see the role of anisotropic parameters. In this paper we 
consider two types of anisotropy and study the effects of 
them on the thermal entanglement. Both the antiferro- 
magnetic and ferromagnetic cases are considered. 

First we briefly review a measure of entanglement, the 
concurrence Q. Let pi2 be the density matrix of a pair 
of qubits 1 and 2. The density matrix can be either pure 
or mixed. The concurrence corresponding to the density 
matrix is defined as 



Ci2 = max {Ai - A2 - A3 - A4, 0} , 



(1) 



where the quantities Ai > A2 > A3 > A4 are the square 
roots of the eigenvalues of the operator 



where the coupling constants J > corresponds to the 
antiferromagnetic case and J < the ferromagnetic case. 
The operators aj± — \ {ujx ± if^jy) {j = 1, 2). The XXZ 
model was initiated by Bethe for the case A = ±1 in 
1931 [g and has been studied for A 7^ ±1 since 1959 [|. 

The eigenvalues and eigenvectors of H are easily ob- 
tained as 
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where l^**) = ^d'^-'^) ^ maximally entangled 

states and |0) (|1)) denotes the ground (excited) state of 
a two-level particle. 

In the standard basis, {|00), |01), |10), |11)} , the den- 
sity matrix p(T) is written as (fc = 1) 
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The square roots of the four eigenvalues of the density 
matrix £112 are 



£'12 = P\2{p\y ® '^2y)p\2^0\y ® '^ly)- 



(2) 



The opeators Ojy{j = 1,2) are the usual Pauli operators 
for the qubit ]. The concurrence C12 = corresponds 
to an unentangled state and C12 = 1 corresponds to a 
maximally entangled state. 

We consider the two-qubit anisotropic XXZ Heisenberg 
model 
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a — (CTi2;Cr2x + criyCr2y + Ao-izCr2z) 

= J((Tl + Cr2- + 0-l-Cr2+) + ^crizO'22, (3) 



Which is the largest eigenvalue depends on the value 
of anisotropic parameter A and sign of J. For antiferro- 
magnetic case ( J > 0) the largest eigenvalue is Ai when 
A < — 1 and A3 when A > — 1. Therefore the concur- 
rences are given by 



1 



Cafm{^) = OforA<-l, 

r< /sinh(^)-e- 
Cafm{^) = max ' 



cosh ; 



for A > -1. 



(7) 



When A — 1, the anisotropic model becomes the 
isotropic model, and Eq.(0) reduces to 



that the concurrences are 1 for different anisotropic pa- 
rameters when T = 0. In these cases the ground state 
is 1 5*^), which is the maximally entangled state and the 
corresponding concurrences are 1. As the temperature 
increases, the concurrence decreases due to the mixing of 
other states with the maximally entangled state. Again 
we see that Tc increases as A increases. 



Ca_fm(1) = max — , 

\e-+3 . 



(8) 
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which is obtained in Ref. [y. From the above equation 
we know that when the temperature is larger than the 
critical temperature Tc = the thermal entanglement 
disappears. For the anisotropic model the critical tem- 
perature Tc is determined by the nonlinear equation 



smh(-) = e T , 



(9) 



For ferromagnetic case (J < 0) the largest eigenvalue 
is A4 when A < 1 and Ai when A > 1. Therefore the 
concurrences are 
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Cfm(A) = for A > 1 
Cfm(A) = m; 
for A < 1. 
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From the above equation we see that no thermal entan- 
glement for the ferromagnetic isotropic Heisenberg model 
(A = 1). The critical temperature is given by the equa- 
tion 



(11) 



From Eqs.(|^) and ( p^ it is find that the thermal en- 
tanglement are same when A = 0. That is to say, the en- 
tanglement exists in the antiferromagnetic and ferromag- 
netic models at the same time. The Heisenberg Hamil- 
tonian with A = is just the quantum XY model. The 
thermal entanglement this model is discussed in a recent 
paper From Eqs.(|^ and ( p^ we also see that the 
concurrences satisfy Cafm{'^) — Cfm{^^)- 

We numerically solved Eqs.(||) and and the re- 
sults are shown in Fig.l. For the antiferromagnetic case 
we observe that the critical temperature Tc increases as 
the anisotropic parameter A increases. Oppositely Tc 
decreases as A increases for the ferromagnetic case. Of 
course the critical temperatures are same when A = 0, 
which corresponding to the XY model. 

Fig. 2 (a) gives a plot of the concurrence as a function 
of temperature for the antiferromagnetic case. It shows 



FIG. 1. The critical temperature as a function of 
anisotropic parameter A for both the antiferromagnetic and 
ferromagnetic cases. The parameter J = 1. 

Another kind of anisotropy is the DM anisotropic an- 
tisymmetric interaction which arises from spin-orbit cou- 
pling 1^,0. Now we consider the Heisenberg model with 
DM interaction 



HdM = ■^[(CTi2;Cr2a; + criyCr2y + AcrizCr22 
+3 ■ (cti X CT2)], 



(12) 



where D is the DM vector coupling. To see the effect 
of the anisotropic parameter D we choose D — Dz and 
A = 0. Then the Hamiltonian H^m becomes 

Hdai ~ —[crix<72x + o-iy<72y + D{aixCr2y — Criya2x)] 
= J[(l + tD)ai+cr2- + (1 - iD)ai-(J2+].{13) 

The eigenvalues and eigenvectors of Hoai are given by 



HdmIoo) = o,HDM\n) = o, 
Hdm\±) = ±jVi + d^\±), 

where |±) = (|01) ± e^®|10)) and 6 = arctanD. 



(14) 
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In the standard basis, the density matrix p(T) is given 
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The square roots of the four eigenvalues of the density 
matrix qyi are 
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We see that the four eigenvalues are independent on the 
angle Q. From the eigenvalues we observe that for both 
antiferromagnetic and ferromagnetic cases the concur- 
rences are given by 
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We see that the entanglement does not depend on the 
sign of the anisotropic parameter D. 
The critical temperature is given by 



Tr 



|J|^/^+J^ 
arcsinh(l) 



1.1346^1 + -0^1 J| 



(18) 



Obviously the critical temperature increases with the in- 
crease of the absolute value of the anisotropic parameter 
£), which can also be seen in Fig. 2(b). The ground state 
of the system is |-|-) or |— ) no matter what the anisotropic 
parameter D is. Fig. 2(b) shows that the concurrences are 
1 when r = 0, which is due to the maximally entangled 
ground state |±). When the temperature is larger than 
the critical temperature the thermal entanglement disap- 
pears. 
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FIG. 2. The concurrences as a function of temperature, (a) 
For anisotropic antiferromagnetic Heisenberg model; (b)For 
Heisenberg model with DM interaction. The parameter 
|J| = 1. 

In conclusion we have studied the effect of two 
kinds of anisotropy on the thermal entanglement in the 
anisotropic XXZ model and the Heisenberg model with 
DM interaction. For the XXZ model it is shown that 
the thermal entanglement exist or not depends on both 
the anisotropic parameters and the sign of exchange con- 
stants J. The thermal entanglement are same for the anti- 
ferromagnetic and ferromagnetic Heisenberg model with 
DM interaction. While in the XXZ model the thermal 
entanglements are different for the antiferromagnetic and 
ferromagnetic cases. In this paper we restrict ourselves to 
the two-qubit case. It is a good challenge to study ther- 
mal entanglement in the multi-qubit anisotropic models. 
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